Quark- meson coupling model with constituent quarks: 
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Abstract 



^ ^ The binding energy of nuclear matter including exchange and pionic effects is 

Q>^ ' calculated in a quark-meson coupling model with massive constituent quarks. 

As in the case with elementary nucleons in QHD, exchange effects are re- 
pulsive. However, the coupling of the mesons directly to the quarks in the 



i nucleons introduces a new effect on the exchange energies that provides an 



extra repulsive contribution to the binding energy. Pionic effects are not small. 
Implications of such effects on observables are discussed. 
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1. Introduction. In a recent series of papers Saito and Thomas (ST) have devel- 
oped a model for nuclear matter and finite nuclei in which the quark structure of the nucleons 
is explicitly considered. The model is based on an original proposal from Guichon P], and 
is known as the quark- meson coupling model (QMC). The model shares important features 
with the Walecka model (QHD) [Q, such as a relativistic saturation mechanism, but the 
incorporation of explicit quark degrees of freedom has nontrivial consequences (for a list 
of references see Ref. 0). Recent developments of the model have been made by Jin and 
Jennings and Blunden and Miller [B] . For earlier studies within the QMC model see Ref. H . 
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Exchange and pionic effects have not been investigated so far in the the QMC model. 
Inclusion of exchange effects (Fock terms) is crucial for assessing the effects of pionic degrees 
of freedom, which are important in view of the connection of the pion to the dynamical 
chiral symmetry breaking of QCD. Since the earlier days of QHD, it has been shown 
that the inclusion of Fock terms requires a sizable renormalization of the parameters of the 
model, as compared to the values fixed at the mean field level. Different parameters of scalar 
and vector mesons imply for instance in a different response of matter to external hadrons 
immersed in medium. Also, the consideration of exchange effects is an initial step towards 
the consideration of short-range correlations (SRC), which have been shown within the tra- 
ditional nuclear many-body theory to play an important role on the saturation properties of 
nuclear matter. 

In this paper exchange and pionic effects are studied in the context of a quark-meson 
coupling model with constituent quarks. The main motivation for using a constituent quark 
model (CQM) is its simplicity. The CQM has provided a good deal of insight into the study 
of a variety of hadronic phenomena at the nuclear scale [Q, and it is natural to expect that it 
should be equally useful for addressing the role of quarks on low-energy properties of nuclear 
matter. In the next section we present the model and discuss approximations. In Section 3 
we discuss the Hartree approximation and in Section 4 we consider the Fock terms and study 
the effects of the pion. Conclusions and future perspectives are presented in Section 5. 

2. The model. The mesons are described by the Lagrangian density of QHD, including 
the 0", uj, TT, and p mesons. We separate from the a, u;, p field operators their mean- 
field, classical values 0O; and po, and obtain a meson quark-quark interaction from the 
fluctuating part of the fields. The exchange of mesons between quarks is treated in a similar 
fashion to the traditional one-gluon interaction [§, with the difference that we do not make 
a nonrelativistic reduction neither of the kinetic energy nor of the meson-quark interaction. 
The exchange energy is the expectation value of the meson quark-quark interaction in the 
state |\E' > of non-overlapping nucleons, which we write in the form of a "Fermi gas of 
composite nucleons" , 





where the nucleon creation operator is given by 
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$p is the three-quark Fock-space amphtude, where p is the nucleon center of mass (cm.) 
momentum, and a the nucleon spin-isospin third components. T^'"^ is a spin-isospin Clebsch- 
Gordon coefficient, a,b,c = (sf) are the quark spin-flavor quantum numbers and eijk is 
the totally antisymmetric tensor in the color indices. The quark creation and annihilation 
operators satisfy canonical anticommutation relations. The largest value of |p| in Eq (|l]) is 
the nucleon Fermi momentum kp, which is related to the nuclear matter density as pn = 
^kp/GTc"^, with 7 = 4(2) for nuclear (neutron) matter. 

Let us briefly consider the nucleon in free space, and collect some useful material for 
later discussions. For simplicity, we neglect gluons and use for the confining potential an 
harmonic oscillator p[, V{r) = — A'^A"/4 With a nonrelativistic kinetic energy, the 
problem is exactly soluble. For the ground state, 

$p(fci, fcs) = Sip -k,-k,- fcs) (^^^ ^ e-'"/' S'<. - ^^■)', (3) 

where b is the r.m.s. radius of the nucleon. In case of a relativistic kinetic energy, the 
problem is not exactly soluble. We then use $ given as above, and treat 6 as a variational 
parameter. In this case, the ground state mass of the nucleon is given by 

3/2 
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where the second term in the square bracket is the energy from the confining potential, 3Cb^. 
The constant C was eliminated making use of the variational condition on M^r. 

In order to calculate the energy density of nuclear matter S we need the interaction 
of the quarks with the unconfined, mean meson fields. We obtain the meson quark-quark 
interaction in a similar fashion to the derivation of the one-gluon interaction in the CQM ^ , 
but with the difference that we use Dirac field operators expanded in a plane-wave basis, 
ij}{x) = /dfc/(27r)3/2 J2sUs{k)e'''-''qif{k), with 



Usik) 



E*{k)+m* I 
\ V E*{k)+m* 



Xs, (5) 



where E*{k) = y^fc^ -|- m*^, and m* = iriq — Qqa (po is determined requiring stability of the 
energy density with respect to variations in 0o- Note that since antiquarks are neglected, 
the three-quark nucleon wavefunction, < 'ilj{xi)'ilj{x2)ip{x3)\B'l{p) >, contains only positive- 
energy components. This ensures that the Dirac Hamiltonian is bounded from below. We 
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note however that neglecting antiquarks is not entirely consistent since it does not lead to a 
well defined Hilbert space, and should be considered as part of the definition of the model. 
The use of a wave function with a complete Dirac structure requires the use of a different 
variational principle to avoid the continuum dissolution problem [|lOl, such as the minimax 



variational principles used in atomic physics for relativistic many-electron systems |TT|]. In 
this case, the energy functional in the space of variational parameters is a minimum with 
respect to the parameters related to the upper component of the Dirac spinor, and a maxi- 
mum with respect to the ones related to the lower component. Such variational calculations 
have been recently applied to quarks models [0, and appear to be particularly suitable for 
nuclear matter calculations like the present one, where a variational calculation seems to 
be imperative. We do not pursue this here, rather reserve it for future elaboration, along 
with the consideration of a proper relativistic structure for the confining potential. This 



last point is still a matter actively studied ||T3[; it is relevant in connection with spin- and 
momentum-dependent effects of the confining interaction, which are not considered in the 
present paper. In the approximation we are using, there is no difference between a potential 
that is the time component of a four-vector and a scalar potential in which all spin and 
momentum-dependent terms are neglected in a nonrelativistic expansion. For a wave func- 
tion with a complete Dirac structure however, a confining vector interaction would lead to 



problems with the Klein paradox, as is well known [Q. One should also note that when Fock 
terms are considered, Eq. is not the most general spinor basis, since exchange effects in- 
duce a momentum dependence for the quark self-energies. In this sense, our approach shares 
similarities with the one of Chin 0] in QHD, where exchange energies are evaluated with 
the mean field self-energies. The same approximation is employed in the Dirac-Brueckner 
approach . 

For nucleons in medium, we use a gaussian form for $p as in Eq. (^), with b replaced by 
a b*, and determine b* variationally. Neglecting effects due to the superposition of the quark 
clusters in medium, such as quark exchanges between different nucleons (these are discussed 
in Section 4) and effects thereof [|l^, the energy density of nuclear matter is given by 



(2vr) 

+ lml4 + lmWo + lmlff,+ ^ Vj^, (6) 
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where is the mean field exchange energy 



4 



+- 



9 f^F dp r^F dp' 



2 7o (27r)3 Jo (27r)3 ^ 



■ujp{p,p)e 



b*'{p'-py/3 



(7) 



where 



'36 



*2 ' 



27r 



Jdkdqe-^b^y^K'^ - + ^) ] A,(t; ^ 

2"(si/i)6ic 



rp(s2f2)b2C2 rp{s3f3)biCi 



V18 V18 Vl8 
where i = p' - p, = (p - t)/3, g = {p' + t)/3, 6E*{k,t) = E*{k) - E*{k + t), and the 
vertices and propagators are 
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where q = — ki = k2 — k^. We have dropped the term proportional to k^k'^ /ml^ ^ in the 
vector potential because of current conservation. 

As the final step, m* and h* are determined through the variational equations d£ /dcl)^ = 
and dS / db* = . Note that 0o is a dynamical variable, and as such can always be determined 
from the thermodynamic argument that an isolated system at fixed baryon number and 
volume (and zero temperature) will minimize its energy. 

3. Hartree Approximation. The Hartree approximation consists in taking = in 
Eq. (BI). The minimization of the energy density with respect to 0o leads to 



3m^ = 3m,-2^7( — 
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Before discussing the numerical results, we start performing a nonrelativistic approxima- 
tion, with the only purpose of getting insight into the problem. All integrals can be done 
analytically. Eq. ([Tl|) becomes 
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and the energy per nucleon (for symmetrical nuclear matter po = 0) 
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The variational condition for b* leads to b*/b = (m^/m*)^/^ > 1, since m*/mg < 1. An 
interesting comparison can be made with QHD, when the Dirac spinors for the nucleons 
are expanded to (9(p^/M^) [|l^. The effective mass of the nucleon, Eq. (11) of Ref. ||T7[| , 
becomes M]^ = Mtv - gl/ml (1 - 3kj,/10M^^ )p^. Comparing this with Eq. above, it 
becomes evident that the effect of internal structure of the nucleon is the term l/2m*^6*^ 



in Eq. (JT^). This term comes from the kinetic energy of the quarks in the nucleon, and is 
of the opposite sign to the QHD term. The internal motion of the quarks provides an extra 
repulsion and acts clS 8b repulsive field in medium. 

We next consider the numerical solutions of the relativistic equations. In this case the 
integrals in Eqs. (|D(0)(|D and (0) must be performed numerically. The determination of 
b* must also be done numerically. We fix the coupling constants such as to obtain a stable 
minimum of E/N — M^v — —15.75 MeV at kp ~ 1.36 fm~^. We use standard values for the 
masses, namely rrig = 350 MeV, rricr = 550 MeV, and = 783 MeV. The parameter set is 
shown in Table I, where we show results for b = 0.6 fm and b = 0.7 fm. 

In Fig. 1 we plot E/N — Mn as a function of kp for different model parameters of Table I. 
In Fig. 1(a), the QHD mean-field result with parameter set (1) in Table I, is the long-dashed 
line. Our Hartree result for b = 0.6 fm, parameter set (2), is the long-short-dashed line. 
Inspection of Fig. 1(a) reveals that our Hartree result provides a nuclear compressibility, 
K, that is smaller than the one obtained in QHD. As seen in Table I, the values for K 
for b = 0.6 fm and b = 0.7 fm are considerable smaller than the value in QHD (MET), 
K = 545 MeV. The qualitative feature of the results is that the effective repulsion provided 
by the internal degrees of freedom of the nucleons increases as b decreases. 

The original QMC model of ST predicted a much smaller mean u field as predicted in 
finite-density QCD sum rules and relativistic nuclear phenomenology |18|. On the other 



hand, when the correction due to the cm. in the ST model is taken to be independent of 
the applied field, the discrepancies are significantly reduced. Our value of g^j is very similar 
to the new value in the ST model 0. 

3. Exchange and pionic contributions. The evaluation of the Fock energy involves 
the summation over the spin-isospin quantum numbers of the quarks and the evaluation of a 
seven- dimensional numerical integration in Eqs. (0)(^). The spin-isospin sum can be easily 
performed using the "substitution rules" |]19[ , which allow to express quark spin-isospin Pauli 



matrices in terms of nucleon spin-isospin Pauli matrices. Regarding the multidimensional 
integral, although it can be carried out using a Monte Carlo integrator, we notice that the 
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integrands of the k and q integrals in Eq. (|^) are concentrated around the values of k and 
respectively . In view of this, we use k for k and q for q in the spinors and in the meson 
propagators and integrate analytically over q and k. What remains is a three-dimensional 
integral, which has to be performed numerically. The minimization of the energy density 
with respect to 0o ^^nd h* is also performed numerically. 

The new values of and g^j are shown in Table I. The readjustment of g^ is relatively 
small, of the order of —10%, and of g^ is of the order of +40%. As in QHD, the effect of the 
Fock terms is repulsive, as seen in Fig. 1(a), where the curve with short-dashes is calculated 
with the Hartree parameter set (2), including the exchange terms from the a and uj mesons. 
However, the internal structure of the nucleon has an extra, nontrivial effect, which can be 
qualitatively understood through a nonrelativistic reduction of the Fock energy from the a 
meson 

(14) 
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The first term in square brackets is the nonrelativistic reduction of the exchange energy 
in QHD [0. The term l/m*^6^ represents the nontrivial effect of the internal structure 
of the nucleon, which is of the opposite sign to the QHD term, and is responsible for the 
enhancement of g^. This term is equivalent to the a dependence of the g^N coupling constant 
in the original QMC work A similar, but smaller effect is seen in the contribution from 
the uj. The enhancement of g^ is of interest for a good description of the spin-orbit splittings 
in finite nuclei p|. 

The vr-quark coupling constant is determined [|T9| by comparing the asymptotic behavior 



of the nucleon-nucleon one-pion-exchange potentials at the quark level and for point nucleons: 
= (9/25)/!,, exp(-m26V2), with fl^ = glMl/^Ml) = 0.98, = 138 MeV and 
gj^^/ATc = 14.4. The effect of the vr is not small, although somewhat smaller than the 
exchange effects from the a and u mesons. The dash-dotted curve is calculated with the H 
parameter set (2), and includes the exchange effects of the a and u. The new readjustment 
of ga and g^^ to saturate nuclear matter at the right place is relatively small, as seen in 
Table I, entry (5). 

Because of the large tensor coupling of the p meson, which is of the comparable magnitude 
to the vr coupling, we also include its effects in the present calculation. The p-quark couplings 
are fixed as in the case of the vr. For comparison with the results of Ref. ||20|| , 



we use 



gj^ /4:7T = 0.55, /np/qnp = 3.7, and nip = 770 MeV. The effect of the p is also repulsive. 



The new values for g^j and are given by the parameter set (6) in Table I. The solid curve 
in Fig. 1(a) is the result of including all mesons. 

As in QHD [Q, the effect of the vr and the tensor coupling of the p are repulsive. This 
is due to the zero-range (ZR) components in the interaction, which would not contribute 
if short-range correlations are taken into account. In order to obtain a more reasonable 
estimate of the effects of vr and p, we follow Ref. |^ and remove this component from the 
interaction. The effect provided by the attractive vr is of the order of 10 MeV, similar do the 
effect in the QHD calculation of Ref. [^]. The long-dashed curve in Fig. 1(b) is calculated 
with the H parameter set (2) and the ZR component of the vr removed, and includes the 
exchange terms of a and uj mesons. The long-short-dashed curve in Fig. 1(b) is calculated 
with the H parameter set (2) and the ZR components from both the vr and the p interactions 
removed. The net attraction provided by the vr and p is of the order of 24 MeV per nucleon.In 
the calculation of Ref. for the same value of fupl QNp = 3.7, the net attraction is of the 
order of 26 MeV per nucleon. 

To conclude, we note that the exchange effects considered here arise from the Pauli 
principle at the nucleon level, i.e. quark-exchange (QE) between different nucleons were 
neglected. QE effects are expected to be small at low densities, but for for nucleons with 
h = 0.6 fm in a medium at the normal nuclear density, they might have interesting conse- 
quences. In particular, they might be of interest for replacing u exchange as the main source 
of the short-range repulsion, in view of possible conceptual difficulties with the coupling of 



an extended meson to quarks inside the nucleons 0]. There is an extensive literature 
in this subject and it is out of the scope of the present paper to review the subject here. 
The present model can be naturally extended to include QE effects following the idea [|T^ 
that the innermost, short-ranged part of the NN force is generated by QE with pion cou- 
phng, and the outer part is generated by meson coupling to the nucleon. In such a picture, 
there is no danger for double counting by taking simultaneously QE and u exchange, since 
the uj coupling to the nucleon is cut-off by the nucleon form factor. We use the approach 
of Ref. |^2[ to derive an effective NN interaction from pion-quark exchange. The resulting 
NN interaction is non-local and gives scattering phase-shifts numerically similar to the ones 
obtained with the resonating group method . After a long calculation, with the extensive 



use of the substitution rules of Ref. |19], the contribution of QE to the energy density of 
symmetrical nuclear matter can be written as 



~?>mlh i2nr Jo {2Tifti ' 



8 



with Ai = 8(3/4)3/2, ^ y^2, A2 = 54, 02 = 0, A3 = 120(12/11)3/2, = 2/33, 
= -44/3, 04 = 1/3, A5 = -272/3(12/11)3/2, as = 8/33. Here, only the contribution 
of the zero-range part of the pion-exchange is shown. The contribution from the long-range 



part is numerically negligible, as also found in Ref. [T^\, where a nonrelativistic form for 
the pion-exchange was used. The new coupling constants (we do not consider p coupling) 
are g"^ = 173.45 and = 56.2. The compressibility is K = 375 MeV and b*/b = 1.07. 
Note that this value of g"^ is very close to the quark-model SU(6) symmetry prediction 
(72/471- ^ 9(^2/477 = 9 X 0.55, which is much smaller than the values used in QHD and in 



one-boson-exchange models []T5|. We have also checked the effect of quark exchange with a 
coupling; the effect amounts to a new readjustment of ga, with no important consequences 
for the values of K and b*/b. 

4. Conclusions and Future Perspectives. The consideration of exchange effects in 
QMC-type of models is the initial step towards the implementation of chiral symmetry and 
short-range correlations. Important consequences are expected when both elements are put 
together, as recently shown by Banerjee and Tjon in a study using density-dependent 
meson masses and couplings. We have shown here that the coupling of mesons directly 
to the quarks has important consequences on the exchange energy, implying in a sizable 
readjustment of the meson-quark coupling constants. A readjustment of coupling constants 
is not meaningless, since different coupling constants imply in different in-medium spin-spin, 
spin-orbit and tensor forces and therefore different responses will be experienced by external 
hadrons immersed in medium. 

An interesting new direction is the replacement of cj-coupling by nucleon overlap at 
short distances. The idea of describing the very short-range part of the NN interaction by 
quark-exchange has long been discussed and we have shown that an interesting consequence 
is a large reduction of the value of g^. Another interesting aspect that is expected to 
influence the saturation mechanism is the inclusion of the low-lying nucleon resonances. 
The present model has the potentiality to incorporate such effects, since with the values 
huj = \lm\b^ = 174 MeV and gqT^ used here it is possible to obtain a reasonably good 
description of the low-mass spectrum of the nucleon, including the Roper resonance |0 . 

An important and necessary next step is to improve on the Lorentz transformation prop- 
erties of QMC-type of models, based on bag or semi-relativistic constituent quark models, a 
problem that is intimately related to the projection of the cm. motion. This is because low 
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energy theorems based on Lorentz invariance impose constraints on the single-nucleon 
energy in matter which have a direct impact on the saturation properties of nuclear matter. 
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TABLE I: Model parameters and results for the compressibility K and b*/b. The parameters 
9u) = '^9qa and = ^Qqa are adjusted to E/A — Mj^ ~ —15.75 MeV at kp ~ 1.36 fm~^. 



Model 


b (fill) 


9 
fJa 





A" (MeV) 


b'-/b 


(1) 


QHD (MFT) 




105.96 


161.10 


545 




(2) 


H 


0.6 


117.00 


90.27 


330 


1.045 


(3) 


H 


0.7 


112.10 


98.80 


340 


1.05 


(4) 


H+EXCH (a + u) 


0.6 


152.89 


88.43 


370 


1.06 


(5) 


H+EXCH {a + uj + n) 


0.6 


165.60 


90.90 


362 


1.07 


(6) 


H+EXCH {a + co + TT + p) 


0.6 


165.66 


89.80 


351 


1.07 


(7) 


H+EXCH (a + + ir^tt) 


0.6 


147.44 


98.75 


353 


1.06 


(8) 


H+EXCH {a + uj + n^tt + Patt) 


0.6 


137.71 


101.50 


340 


1.05 




FIGURE 1. The total energy per particle as a function of the Fermi momentum, (a) The 
long-dashed curve corresponds to QHD (MFT). The long-short-dashed curve is our Hartree 

result, parameter set (2) in Table 1, and the solid curve is for parameter set (6). The short- 
dashed (dot-dashed) curve is calculated with the H couplings (2) including exchange terms 
of cr + a; (cr + a; + tt). (b) The solid curve is for parameter set (8), with the zero-range parts 
of TT and p removed. The long-dashed (dash-dotted) is calculated with the H couplings (2) 
and exchange terms of o" + a; + TTatt {cr + cu + TTatt + Patt)- All curves are for b — 0.6 fm. 



12 



